Transport in vortex state of d-wave superconductors at zero temperature: 

Wiedemann-Pranz violation 
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We show that the Wiedemann-Franz law is violated at zero temperature in the vortex state of 
a d-wave BCS superconductor with isotropic impurity scattering. We use a semiclassical approach 
to include the Doppler shift experienced by the quasiparticles due to the circulating supercurrents 
and consider as well the Andreev scattering from an array of vortices assumed to be randomly 
distributed. We also show that the vertex corrections to the electric conductivity, which can be 
large when there is significant anisotropy in the impurity scattering, become unimportant as the 
magnetic field is increased. For the thermal conductivity, the corrections remain negligible as in the 
absence of a magnetic field. 

PACS numbers: 74.25.Fy, 74.72. Bk, 74.25.0p 



I. INTRODUCTION 

The Wiedemann-Franz (WF) law is a hallmark of 
metallic behavior. This law states that the ratio of the 
thermal conductivity (k) to the dc optical conductivity 
((t), is given by a universal number times the tempera- 
ture (T). The Lorenz number is defined by L = n/{Ta), 
and is expected to equal 7r^/(3e^) (we use h = kB ^ !)• 
This law is well obeyed in most metals, and can be theo- 
retically justified well beyond the case of a simple metal. ^ 

Recently several experiments have been designed to 
test the validity of the WF law, particularly at very low 
temperature, in the cuprate superconductors. The law is 
well obeyed in overdoped Tl2Ba2Cu06+5 which is forced 
into its normal state through application of a magnetic 
field. ^ A similar result was also found in optimally doped 
Bi2+xSr2_a;Cu06+5i^ These results clearly favor the ex- 
istence of ordinary quasiparticles in these materials and 
do not support the spin-charge separation picture of the 
Luttinger liquid. "^'^'^ In this picture collective modes re- 
place the quasiparticles with spinous carrying the spins 
and holons the charges, and the WF law does not holdii 
On the other hand the normal state of a slightly under- 
doped sample of Pr2_a;Cea;Cu04_y studied by Hill et al& 
revealed serious violations of the WF law, which indicates 
the existence of an exotic normal state. The data puts 
severe constraints on the state involved. For example, 
formation of a competing hidden ordered state such as 
the d-density wave (DDW) state is not expected to lead 
to a violation of the WF lawiS 

In this paper we will confine our discussion to the su- 
perconducting state with a d-wave gap symmetry. In par- 
ticular we will study the mixed state where vortices are 
created by the application of an external magnetic field 
H perpendicular to the two dimensional Cu02 plane. 

Near the critical temperature Tc, the cuprates exhibit 
strong inelastic scattering. The inelastic scattering rate 
is of order a few times Tc.^^ The large peak seen in the 
microwave conductivitjiiiii^ as a function of temperature 
is widely interpreted as evidence for the collapse of the in- 



elastic scatteringiiii as T is reduced below Tc- The idea 
is that the development of a superconducting gap sup- 
presses the low frequency density of bosons which causes 
the superconductivity. This suppression is a general char- 
acteristic of any electronic mechanism, in which the bo- 
son is a collective mode of the superconducting electrons. 
Consequently for T <C Tc elastic impurity scattering will 
dominate even in the purest samples in which the quasi- 
particle mean free path can be of order one micron^ 
Therefore in the low temperature regime, which is the 
region of interest here, it is sufficient, as a first approxi- 
mation, to ignore the specifics of the mechanism involved 
in the inelastic scattering. A d-wave BCS formulation 
which includes impurity scattering will suffice. 

For a d-wave gap symmetry both the quasiparticle 
density of states and the effective impurity scattering 
rate acquire an important frequency (w) dependence. 
These essential ui dependences can lead to new phe- 
nomena. For example impurities modify, in an essen- 
tial way, the quasiparticle density of states N{uj) in the 
limit u! 0-^^'^^ In the pure case the density of states 
approaches zero linearly; however, with impurity scat- 
tering it takes on a finite value at zero frequency which 
leads directly to the concept of the universal value of the 
conductivity at u; = and T — 0, independent of im- 
purity content fii This behavior also holds in the case of 
the thermal conductivity and has been verified in many 
experimentsiiSii^iSSiSiiS^ However, the WF law still holds 
in this limit. 

If finite temperature is considered, the WF law is 
violatedfSS Violations can be either positive or negative 
(the Lorenz number is greater or less than one), depen- 
dent on the details of the impurity potential. In the uni- 
tary limit, the Lorenz number can increase rapidly with 
increasing T while in the Born limit (weak scattering) it 
can decrease. In pure samples it drops to one half its 
original (T — 0) value within a few degrees Kelviui^ In 
the simplest approximation, k and a both depend on the 
overlap of thermal factors with a product of the density of 
quasiparticle states N{uj) and the scattering time t{llj). 
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The thermal factors for the electrical case are peaked 
about oj = while for the thermal case they peak at 
higher frequencies of order T. This means that a different 
frequency range of N(lu)t(lu) is most importantly sam- 
pled in the two cases and the WF law no longer holds at 
finite temperaturcf^i The Lorenz number can be smaller 
or greater than its conventional value depending on the 
behavior of N{uj)t{uj) as a function of u). For the su- 
perconducting state this important function N{uj)t{u!) is 
different for the electrical and thermal conductivity be- 
cause of differences in coherence factors but the physics 
described above for the normal state continues to hold. 

Recent extensions of the theory at T = to include 
vertex corrections, which is necessary when the impu- 
rity potential is anisotropic, have led to modifications of 
the above picture. Durst and Lee (DL) found that in 
these circumstances there can be important corrections 
to the electrical conductivity while the thermal conduc- 
tivity is almost unaltered.^''' Consequently the WF law 
no longer holds even at T = and the Lorenz number is 
found to be reduced. The amount of reduction depends 
on the impurity potential Vi (Ref.^^) and is largest for 
the unitary limit {Vi oo). The physics underlying this 
phenomenon is well described in Ref.^^ and has its ori- 
gin in the fact that the thermal conductivity depends on 
the group velocity of the quasiparticles while its electrical 
counterpart is sensitive to their Fermi velocity. 

In this paper we want to consider the effect that a mag- 
netic field (H) applied perpendicular to the two dimen- 
sional d-wave superconducting plane will have on the WF 
law and on the vertex corrections for the in-plane trans- 
port. It is well known2L22i224^ that the application of H 
will increase the density of states at zero frequency iV(0) 
and thus increase both a and k. When H is larger than 
the lower critical field Hd but still much less than the up- 
per critical field Hc2 , vortices form and these will provide 
supercurrents as well as an added scattering mechanism 
for the quasiparticles. Here we treat effects of the super- 
currents semiclassically and this additional scattering in 
a phenomenological model as Andreev scatteringSL^Si^^ 
which can be different for charge and heat transport. We 
found the WF law is violated at T = in the mixed 
state of a BCS d-wave superconductor. This result holds 
even when the vertex corrections for anisotropy are not 
taken into account, i.e. at the level of the bare bubble 
diagrams in the treatment of the correlation functions 
involved. We have also studied the influence of an exter- 
nal magnetic field on the vertex corrections. For the case 
of the electrical conductivity, it is found that the vertex 
corrections present at zero H can be rapidly suppressed 
as H increases. For the thermal conductivity, the vertex 
corrections are never important. 

This paper is organized as follows: In Section II, we 
describe the formalism for d-wave superconductors with 
impurities in the vortex state. As a first approximation, 
Matthiessen's rule is applied to consider the impurity and 
vortex scattering. We also study self-consistently, in Sec- 
tion III, effects of the magnetic field on the impurity scat- 



tering, which is beyond Matthiessen's rule. In Section IV, 
magnetic field effects on the vertex corrections to a and 
K is presented. Section V is devoted to conclusions. 

II. THEORETICAL APPROACH 

The calculation of the thermal and dc conductivity 
proceeds in the standard way through the evaluation of 
the appropriate current-current correlation functions for 
heat and charge current, respectively. Formally, the re- 
quired correlation function is written in Nambu notation 
(2x2 matrix notation) and an analytic continuation from 
the Matsubara to real frequencies is needed before the 
zero frequency limit is taken. The superconducting ma- 
trix Green's function in the Nambu notation takes the 
form G(k, iLUn) = {^to + A^n + ^kT3)/{{:;^ - - ), 
where f's are the Pauli matrices in spin space, Ak a d- 
wave order parameter, and is the electronic energy 
dispersion in the normal state. The correlation function 
can be expanded into a product of two Green's functions 
and a vertex function. At the lowest level of approxima- 
tion, the vertex is neglected and only the bare bubble is 
retained. While we start with this approximation later 
we will generalize our work to include the vertex correc- 
tions as in RefiSi; these vanish when the impurity scat- 
tering is isotropic, but are required when the scattering 
is anisotropic. 

To include the effects of an external magnetic field H, 
for Hci ^ H Hc2, we employ the semiclassical ap- 
proximation which includes, in the single particle Green's 
functions, the Doppler shiffe2L22i22i^2iM associated with 
the circulating supercurrent around a vortex core. This 
is accomplished formally by changing the Matsubara fre- 
quencies iLUn in the single particle Green's function to 
iojn — Vs • k, where is the velocity field of the super- 
current and k is the momentum of a quasiparticle on the 
Fermi surface. Now the new Green's function depending 
on Vs • k is inhomogeneous in space variable r through 
the dependence of the velocity field Vs(r) on position 
r measured from the core of a vortex. Thus thermal 
and transport coefficients become local quantities which 
need to be averaged over a vortex unit cell, assumed 
to be a circle of radius R for simplicity. The radius R 
is related to the coherence length and Hc2 through 
2R = ^o\/27ra~^ {Hc2/ H)'^^^ , where a is a geometrical 
factor of order unity. The average of any quantity J-{r) 
is then 

^(H) = ^ / Hr) • (1) 

We will assume that S^o R and that the circulating cur- 
rent velocity field can be taken to drop as 1 /r throughout 
the vortex unit cell. For a magnetic field applied perpen- 
dicular to the Cu02 plane, the Doppler shift Vg ■ k takes 
the form Eh/ p with p — r/R and av f /2yJi:H /^^q, 

where Vf is, the Fermi velocity and <f>o is the flux quantum 
related to Hc2 by $o = 27r^0'^c2- 
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Impurities can be introduced into the calculations 
through the impurity self energy term which is to be eval- 
uated self-consistently in a t-matrix approximation. The 
retarded self energy S^^ret is given by 



(2) 



data as done in our previous worki^ This linear depen- 
dence of the vortex scattering on Eh implies that the 
vortex cores behave like scattering centers. This means 
that we consider two effects due to a vortex; the Doppler 
shift due to the supercurrent as well as the vortex scat- 
tering due to the core. The total scattering rate denoted 
by 7tot(t^) is then given by 



where Go{l^) = [2^7V(0)]-i Tr[G,e*(k, c^)] with Gret 
the retarded Green's function. The summation over k 
is to be carried out over the two dimensional Brillouin 
zone. In Eq. |21 T is related to the impurity density rii 
by r = rii/ [7rA^(0)] and c is the strength of the impurity 
potential Vi] specifically c — l/{nN{0)Vi). The unitary 
limit corresponds to Vi —>■ oo, i.e. c — > 0, while the 
Born limit corresponds to Vi — > (c — s- oo). As we have 
done in our previous work on the thermal conductivity,'^^ 
in the mixed state we assumed as a first approximation 
that Si, ret remains unmodified by an external magnetic 
field. That is, appealing to Matthiessen's rule, we add 
on to the term ImSi^^et = li^) a term which describes 
phenomenologically the Andreev scattering of the quasi- 
particles off the vortex core. A simple expression for the 
quasiparticle-vortex scattering has been given by Yu et 
alm^ Examination of their formula for the geometry used 
here {H _LCu02 plane) reveals that it is proportional to 
the magnetic energy Eh as also argued through dimen- 
sional analysis in Refi^ The constant of proportionality 
b is to be determined by comparison with experimental 



7tot(t^) = 7(^) + bEn , 



(3) 



where the first term 7(0;) is due to impurities and the 
second term bEn is due to vortex scattering. Later we 
will include the effect of H on the impurity scattering 
itself. 

A simplification in the required algebra results when 
the limit of low T is considered. For T <C Aq (am- 
plitude of a d-wave gap), it is appropriate to make a 
nodal approximation^^ for the summation over the Bril- 
louin zone. Application of this approximation to Goi^j) 
in Eq. © leads to 



GoiCo) 



Ao 4Ao 



(4) 



where a cut-off of 4Ao has been applied to the real part 
of Go in Eq. In a similar way, application of the 

nodal approximation to the expression for the thermal 
conductivity k{T,H) leads to 



k{T, H) 
T 



(5) 



r 



where ^/(g) is the Fermi (gap) velocity and /(cjj) = 
[1 -I- exp(a;/T)] ^ is the Fermi-Dirac thermal factor. The 
thermal conductivity in the plane = k/2 in the nodal 
approximation. This relation also holds for the electrical 
conductivity. In deriving Eq. Q we have introduced a 
vortex distribution function V{f) through the equation 



V{e) 



1 



dr(5(e-v,(r)-k) 



(6) 



where the integration is over the vortex unit cell. The 
actual thermal conductivity is obtained by averaging the 
local quantity over the unit cell. In Eq. Q, the function 
Ak{u}, e) can be written as 



A^{Lo,e) = 1 + 



7t('^) 



arctan 



(7) 



where Q — uj — Re [Sj,et]- Several choices can be made 
for the vortex distribution function and results depend 
somewhat on the choice. However, qualitative results are 



not sensitive to the choice. To be specific we report here 
the results based on the Gaussian distribution 



V{e) 



1 



■ exp 



R2 



(8) 



Note that for = we can take e = in Eq. O , namely, 
A^{oj, 0). If we further take T ^ 0, we obtain 

2 



^ l IVf ^ Vg 
T 3 \Vg Vf 



(9) 



This value is the well-known expression for the universal 
limit of the thermal conductivity independent of the im- 
purity scattering rate at zero frequency (700)- If instead 
we had taken T first in Eq. jSJ), we would have 



lim 



«(r,g) 

T 



= + I deV{e)AMe) , (10) 



J. 



where 
A(0,e) = 1 



7t(0) 



7t(0) 



arctan 



7t{0), 



(11) 
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Note that in this case only the zero frequency hmit of 
7( enters. At any finite T, however, the co dependence 
of 7t(w) becomes very important, and the T dominated 
region T ^ (0) can be quite different from the impurity 
dominated region for which T <C 7t(0). FinaUy we note 



that taking H ^ in Eq. H10(l again gives the value of 
the universal limit as it should. 

Similar algebra gives the expression for the zero fre- 
quency value of the optical conductivity <j{T, H), which 
takes the form 



(12) 



r 



where 



1 



u! — e I LO — e 

■ arctan 



7t(^) 



(13) 



In Eq. H12|) . e is the electron charge. We note that the 
difference between An and Aa is traced to the different 
coherence factors that enter, respectively, the electrical 
and thermal conductivity. In terms of the spectral func- 
tions v4(k, w) and _B(k, lj) associated, respectively, with 
the diagonal and off-diagonal part of Green's function, 
the optical (thermal) case involves yl^(k, lj) ± i3^(k, cj), 
respectively iS As T — > 0, we obtain 



(14) 



If we also take = 0, we recover the well-known result 
of the universal value of the dc conductivity 



cjoo = — 



(15) 



Now for H — and T ^ 0, the Lorenz number L = n/Ta 
can be evaluated to test the WF law. Using Eqs. © and 
(|15|l . we obtain ioo = kqo/ (Taoa) ~ tt^/ (3e^) assum- 
ing Vg <C Vf. This is the universal constant normally 
associated with the WF law. 

Remaining at T = we can consider the effects of 
a magnetic field H on the WF law at the level of bare 
bubble i.e. ignoring the vertex corrections. The normal- 
ized thermal (k/koo) and dc (cr/croo) conductivity can be 
written as 



K 

and 



1 



dx e 



1 



X 



arctan 



Ft 



= —= I dx e 

o^oo V' 



X fx 
1 -I- — arctan — 
It V 1 1 



(16) 
(17) 



where Fj = 74 (0) / Eh = 700 /Eh + b with b = for the 
thermal conductivity and b ^ ba for the dc conductivity. 
These two quantities are not the same for Andreev scat- 
tering off a vortex core. On physical grounds we expect 



bn > ba because Andreev scattering is more effective for 
heat transport than for charge transport. Evaluation of 
these parameters requires a microscopic mechanism for 
the quasiparticle-vortex scattering, which is beyond the 
scope of this work. Here we effectively treat them as phe- 
nomenological parameters and present results for various 
values of b. The value of 6^ has been determined for two 
samples in our previous work2& on the thermal conduc- 
tivity. In the ultra pure sample of Hill et al.^ its value is 
about 6k = 0.3 while in the impure sample discussed by 
these same authors we found that Andreev scattering is 
not important compared with the impurity scattering so 
that the value of b^ is effectively zero. 

A first result can be obtained analytically from 
Eqs. H16(l and (|17(l assuming no vortex scattering {b^ = 
6a = 0). In the high field limit Eh > 700 or F* < 1 
we can expand k/kqq and ct/ctoo in terms oi "/qq/ Eh and 
obtain 

— = IF- V(21nrt-l-7e-l)rt (18) 
kqq 41 t 4 

where 7e — 0.577 is the Euler number, and 



cr 

Coo 



2Tt 



(19) 



Note that Eh ^ y H indicating that k/koo ^ VH + 
O {\n{H)/^/H^ while ct/ctoo ^ ^/H + O (l/Vi?) . The 
Lorenz number is given by (Ft << 1): 



(20) 



where Lqq is the universal value of the Lorenz number. 
This means that the Lorenz number at a high field is half 
its universal value, which is a serious violation of the WF 
law at zero temperature. We should emphasize that to our 
knowledge this law has always been obeyed in the past at 
T — 0, even in an exotic state such as the DDW stateA^ 
Our result represents a significant exception. 

For finite but small vortex scattering (6^ ^ 1 and b^ ^ 
1), k/koo and a/aoo both saturate to values determined 
only by 6„ and 60-, respectively, in the limit Eh » "foo, 



K 
KOO 



46, 



[21n(6,)-|-7e-l] 6, 



(21) 
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Fig.l ( Kim et al ) 
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FIG. 1: The Lorenz number L normalized to its universal 
value Loo — ■'r^/(3e^) as a function of a magnetic field H 
applied perpendicular to the Cu02 plane. The top frame is 
for a pure sample (BZO) with 700 = 0.3/s' while the bottom 
frame is for a dirtier sample (YSZ) with 700 = I'iK. (See 
Refsi— *— for details.) Various choices of the vortex scattering 
parameters 6„ and ha are considered: {bK,ba) — (0.3,0.3) for 
the solid curves in both frames, (0,0), (0.3,0.1), and (0.3,0) 
for the dashed, dot-dashed, and dotted curves, respectively. 



and 



As bf; and b„ 



a 

CTOO 



26. 4 
0, we obtain 



L 



26k 



(23) 



Since we expect 60- to be less than the Lorenz number 
should be smaller than its universal value in this regime 
and Eq. 123|l could be used to find the ratio ba/b^ of 
the two parameters associated with the vortex scattering. 
Note that for 6„ = 60- — > 0, we recover the previous result 
of L/Lqo = 1/2. In our previous work^S, we were able to 
determine b^ from the saturated value of the thermal 
conductivity as a function of a magnetic field, and so 
Eq. H23I) would give 5. directly from the Lorenz number. 

We performed numerical calculations of the Lorenz 
number for various values of 6^ and 60- at T = 0. In 
Fig. 1, we plot results for L/Lqq as a function of H in 
two specific cases. The top frame is for 700 — O.Sif, a 
value of the impurity scattering rate deduced from con- 
sideration of the thermal conductivity of an ultra pure 
sample35-36 of YBCOe.gg denoted by BZO and the bot- 
tom frame is for a dirty sample denoted by YSZ with 



700 = 13K. The values of {bii,bc) are shown as labels 
to the curves. For example, (b^^b^) — (0.3,0.3) for the 
solid curve. For the dashed, dot-dashed and dotted curve 
{b^,ba) = (0,0), (0.3, 0.1), and (0.3, 0), respectively. Note 
that the analytic result Eq. H20() is verified by the dashed 
curve with (0,0) in the top frame. In fact, in this case, 
the value 1/2 is reached already around H = 1 Tesla. 
For the bottom frame (YSZ sample) the reduction to 
1/2 is not yet seen even at 10 Tesla because the satu- 
ration is controlled not by H but by the ratio Eh/jqq, 
which is needed to be much larger than 1 in order to 
realize Eq. 1^0)1 . Estimates of En/^foo are about SO^H 
(Tesla" ^) for the top frame which is much larger than 1 
for = 1 Tesla. But the bottom sample is dirtier by 
a factor of approximately 40. Eq. (|23|l is also verified in 
our numerical work. Next we turn to the consideration 
of the effects of a magnetic field on the impurity scatter- 
ing. This takes us beyond the simple Matthiessen's rule 
which we relied on in this section. 



III. MAGNETIC FIELD EFFECT ON 
IMPURITY SCATTERING 



We now return to the impurity self energy given by 
Eq. 10). So far, for simplicity, we have treated the self en- 
ergy approximately, by appealing to Matthiessen's rule; 
we simply added the Andreev scattering effect bEn to 
the impurity scattering 7 which is itself left unaltered 
(22) by the magnetic field. Here b = 6k(ct) for the thermal 
(dc) conductivity as we mentioned. Of course, in reality 
the existence of a vortex lattice, which introduces spatial 
inhomogeneities, will modify the impurity term itself in 
Eq. We need to return to Eq. |(2Jl for the retarded 
self energy and include the effects of the magnetic field 
on Go itself. In obtaining Eq. (j^J an impurity average 
has been performed in the standard way to restore trans- 
lational invariance. After this procedure is carried out, 
the momentum becomes a good quantum number again 
and one can use the semiclassical approximation to take 
into account the Doppler shift. This re- introduces spatial 
inhomogeneities. A second average over a random vortex 
distribution is then taken. The legitimacy behind such 
separate averaging procedures is an assumption that the 
impurities and vortices are uncorrelated. 

To deal with both impurity and vortex averaging in 
a manageable way, what we propose here is to replace 
Go(w) in Eq. (j^J with its vortex averaged counterpart. 
Returning to Eq.Q valid at low T, let us introduce the 
Doppler shift and take an average over, for example, the 
Gaussian distribution. Then in the Born limit, we ob- 
tain for the total scattering rate 7tot(<j-') at a; = using 

^t,ret{EH) = -l-liEn) 
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Itot = 477B / deV{t) 



Itot In 



e arctan 



hE 



, 7tot , 



(24) 



where 775 = n^V^ / {2TTVfVg). The factor 4 in Eq. 
appears for isotropic impurity scattering. When we 
consider the anisotropic scattering, this factor wiU be 
changed to reflect the scattering anisotropy. Eq. H24|l is to 
be solved self-consistently for the effective total scattering 
rate jtotiEn) = i{Eh) + bEn- Results for the impurity 
scattering ^{Eh) are plotted in Fig. 2, where ^[Eh) is 
normalized to the zero frequency and zero field scatter- 



ing rate 7(0,0) = 700- The magnetic energy in Fig. 2 
is also normalized to 700- The top frame applies to the 
Born limit and the bottom to the unitary limit. Three 
cases are considered. The solid curve is for 775 ~ 0.04 
and b = 0.1, the dashed curve for rjB = 0.08 and b ~ 0.1. 
The dotted curve is for & = and applies to all values 
of rjB , as is shown below. After some simple algebra we 
find from Eq. if^ 



477s / deP{e) 



Itot In 



e arctan 



, Itot , 



bEn = 



(25) 



where e = e/700, 7*04 = jtot/ loo, and Eh = EhMoo- It 
is clear from Eq. 1)25(1 that for 6 = 0, the resulting equa- 
tion does not depend on the value rjB- Consequently, 
the dotted curve in the top frame of Fig. 2 results re- 
gardless of a specific value of 775. For any finite b, how- 
ever, the result does depend on rjB- We note that in 
all cases in the Born limit, the effect of the magnetic 
field H is to increase the scattering rate, and that for 
a given value of Eh /joo the increase is largest for the 
smallest rjB value. These results confirm the expecta- 
tion based on the following argument. Let r^r be the 
impurity scattering rate in the normal state. For the 
Born limit, the zero frequency scattering rate in the su- 
perconducting state is much less than Tn while for the 
unitary limit it is larger. In fact, for the Born limit the 
decrease in 700 is exponentially dependent on Fjv • As we 
introduce more vortices by increasing the magnetic field, 
more normal regions appear at the expense of supercon- 



ducting ones and, hence, we would expect that, on aver- 
age, the effective scattering rate would increase towards 
Fat. This is precisely what we have found. On the other 
hand, in the unitary limit 700 is very much larger than 
Fat. a useful and well-known approximate formulai^*^ 
for a relation between these two quantities in this limit is 
7oo — O.GSvFAq- Thus, as H is increased we expect now 
that J (Eh) will decrease because it should tend towards 
its normal state value. This expectation is confirmed in 
the bottom frame of Fig. 2, where we present numerical 
results for two values of b; namely, b ~ and b = 0.1, 
and two values of rjij — 2TTniVfVg/'^QQ (12 and 20). One 
can estimate rju using F = ni/TTN{0), 700 — O.GSVTAo, 
and Aq ~ T:N{0)vfVg in the nodal approximation. In all 
cases, as expected ^{Eh) decreases and the decrease is 
larger when 6 7^ 0. The results in the unitary limit are 
based on the following equation: 



Itot 



Itot In 



Vu 



po 



bE 



H 



e arctan 



Itot 



(26) 



One can see that, in this case, even for 6 = 0, ^{Eh) still 
depends on rjij. These expressions are used to obtain 
results at the bare bubble level for the dc and thermal 
conductivity as a function of Eh (see Figs. 3 and 4 below 
for 6 = and b = 0.3, respectively). 

Even when Andreev scattering from randomly dis- 



tributed vortices is not considered, we still see, in Fig. 3, 
that including the effect of H in ^tot increases both k 
(top frame) and a (middle frame) over the values (dot- 
ted curve) it would have if 700 was used instead of j{Eh) 
(i.e. no effect of H is included in 7) in the unitary limit 
while the opposite holds in the Born limit. This is as ex- 



Fig. 2 ( Kim et al ) 



Fig. 3 ( Kim et al ) 
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2 4 6 8 



FIG. 2: The magnetic field dependence of the impurity 
scattering rate j{Eh) as a function of the magnetic energy 
Eh oc ^/H. The normalization factor for 7 and Eh is 700 • The 
top frame is for the Born limit with rjB = riiV,^ /{2-irVfVg). The 
values of 77s are 0.04 and 0.08 as labelled. The vortex scatter- 
ing parameter is chosen to be or 0.1. Note that for 6 = 0, the 
value of 77B is irrelevant for -^((Eh) in the Born limit. The bot- 
tom frame is for the unitary limit with jyr/ — li^niVfVgjy^-^Q. 
The values of r\u are chosen to be 12 or 20, and 6 = or 0.1. 



pected from Fig. 2, where it was shown that H reduces 
^(Eii) in the unitary Umit while it increases '-/{Eh) in the 
Born limit. Further, in Fig. 4, where a non-zero value of 
b = 0.3 is considered, the effects remain as in Fig. 3 but 
all variations with Eh are much reduced. This means 
that Andreev scattering off vortices offsets to a large de- 
gree the increased conductivity due to the increase of 
quasiparticles created by the Doppler shift. Returning 
to the bottom frames of Figs. 3 and 4, we note that in 
all cases the Lorenz number is reduced below its univer- 
sal value Loo as the field is increased. In particular for 
6 0, the decrease is very rapid, and at Eh /"/qo = 2 the 
unitary case has dropped by more than 30% while in the 
Born limit it is by more than 20%. Also for this ratio, 
including the effect of H on the impurity scattering is not 
important in the unitary limit but is more important in 
the Born case. This also holds for b = 0.3 shown in the 
bottom frame of Fig. 4. There is no qualitative difference 
between the 6 = and b = 0.3 cases but the variations 
are not as rapid when Andreev scattering is included in 
the calculations. 



IV. VERTEX CORRECTIONS TO 
CONDUCTIVITY 

We now consider the effect of vertex corrections to the 
dc and thermal conductivity of a rf-wave superconductor 

in the vortex state. This problem was treated in de- 
tail by DL in the case when the impurity scattering is 



b b 

- iL^lz" 

- Tie s 0.05 


1 


1 , , . 
Unitary ^^^^ 






Born 






Unitary . 






Born 


Unitary"" 
1 


1 


Dotted curve if 7 (E^) = Yqo _ 



1.5 1 — • — I — • — I — I I ' I 

2 4 6 8 10 



FIG. 3: The thermal (top), dc (middle) conductivity and 
Lorenz number (bottom frame) as a function of the magnetic 
energy Eh- All calculations are based on the bare bubble di- 
agram; however, the magnetic field dependence of the impu- 
rity scattering is included. The vortex scattering is neglected 
(6k = ha =0). The solid curves in all frames are for the uni- 
tary limit (r/f/ ~ 12) while the dashed curves are for the Born 
limit (r/B = 0.05). The dotted curves arc for the case when 
the magnetic field dependence of the impurity scattering is 
neglected; namely, 700 is used instead of 'yiEn) in the self 
energy. 



Fig. 4 ( Kim et al ) 




FIG. 4: The thermal (top), dc (middle) conductivity and the 
Lorenz number (bottom frame) as a function of the magnetic 

energy at the bare bubble level including the field dependence 
of the impurity scattering. In this figure, the Andreev scat- 
tering is considered (6^ = 60- = 0.3). Other parameters are 
the same as in Fig. 3. 
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anisotropic but no magnetic field is applied to the sam- 
ple. The necessary mathematics is rather lengthy so that 
here we will deal only with the essential modifications to 
the formalism of Ref.^^ that are required when a mag- 
netic field is applied. Following Durst and Lee, we also 
introduce impurity scattering potentials to the same node 
(Vi), the adjacent node (V2), and the opposite node (V3). 
Thus the anisotropy can be parametrized by R2 = V2/V1 
and i?3 = V3/V1. If i?2 = ^3 = 1, the impurity scat- 
tering is isotropic. Our aim here will not be to give the 
most general formulation but rather consider a simpli- 
fied case which can nevertheless provide a preliminary 
understanding of the effects of the magnetic field on the 
vertex corrections to the conductivity. The most impor- 
tant modification in the formalism is to replace, for all 
relevant calculations, the Green's function by its aver- 



age over the vortex distribution. This replacement may 
limit the generality of our results but allows us to get a 
concrete numerical evaluation of the resulting vertex cor- 
rections. A central function in the calculations is F{z) 
as in Refi2^ where 



F{z) = -AnvfVg 



i^nyvfVg 



G{p, z) (27) 



with z = iLO — Yii — Suor- Note that Suor is the self 
energy due to vortices so that it includes both effects of 
the Doppler shift and Andreev scattering. 

Let us consider the vertex corrections to the dc con- 
ductivity. In terms of F[z) the dc conductivity can be 
written as 



J 



avc 1 



Coo 



^ i 1 deV{e) Re 



io). 



I- A 



(0) 



I- A 



(0) 



(28) 



r 



which replaces Eq. (3.13) of DL. In the above Eq. (^5)) . 

(29) 



oz 



and 



F(z)-F{z*) 



with zo = — e -l- i7t(0) and 
rPa 



F{z) 



pdp- 



2z m — 

z 



(30) 



(31) 



The parameters A^^ and ^2"^ are due to the scattering 
anisotropy. They will be specified later but for isotropic 
scattering they vanish. In this instance, we obtain 



1 



— = - / deV{e)Rc l!f\e)-ir{e) 

CTOO 2 



(32) 



which needs to reduce to the expression we obtained in 
the previous section where the vertex corrections were 
neglected. To see if we recover the bare bubble results 
we need to calculate and 



/(")(,) ^21n^- 2 
Zo 



and 



2 / ipa 
— — Im <^ zq In 

7* L zo 



(33) 



(34) 



Now we obtain 
where 

arg 



7t(0) 



2 — 2 — TTT-arg — — 2i arg — 



Zo 



Wo 



arctan 



Zo J V7t(0) 
Thus the dc conductivity becomes 



a 

1700 



deVie) 



1 



7t(0) 



■ arctan 



7t(0) 



Zo 
(35) 



(36) 



(37) 



which is what we had before at the bare bubble level 
(Eqs. (O and 113). 

Now we consider the vertex corrections to the dc con- 
ductivity in the Born limit. In this case, following DL 



we obtain after the appropriate generalization A\ 



(0) 



and <^ = a'^'ilD with a 



1(0) 



(o)/r(o)\ 



v(o) - JO) 



{O^JdeV{e)ll''>{e) (z = 1,2) and 

r2 



Note in this case, 
Po 



' (l-i?3) 



AirvfVg 



In 



7(0) r,B{l + 2Rl + Rl) 



where 



(38) 



(39) 



where (l + 2R^ + i?|) 4 for the isotropic scattering 
(i?2 = -R3 = 1) as we mentioned for Eq. (|24|) in the 
previous section. The expression for the dc conductivity. 
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including the vertex corrections aye, can be simplified 
and becomes 



Coo 



{if) - {ID 



(0)\ 



Since 



(0)n 



we finally obtain 



(40) 



1 



1 



(41) 



For the numerical calculations we choose Ri = 0.9 and 
i?3 = 0.8. Since it is assumed that V\ is small, we retain 
only terms of order . In Fig. 5, we show numerical 
results with r\B — niV^ /2'KVfVg — 0.05. In a rigorous 
sense i^b since Vi — > in the Born limit; however, 
this value is unrealistic for an actual sample. In the top 
frame of Fig. 5, we plot crvc/o'oo as a function of Eh/ Joo 
for bcr — (solid curve) and ba = 0.1 (dashed curve). At 
small values of Eh , the Andreev scattering does not have 
a large effect on cyc/o'oo- For large values of Eh, how- 
ever, the difference can be much larger. The conductivity 
also rises significantly above its universal value. This rise 
is not to be assigned to the vertex corrections. In the bot- 
tom frame we plot uyc'/cr, where a is the dc conductivity 
at the bare bubble level. We see now that in both cases, 
including {b^ ^ 0) and ignoring {b„ = 0) Andreev scat- 
tering, the application of a magnetic field decreases the 
effect of the vertex corrections although the decrease is 
not large in magnitude. Note that ior H — one can get 



analytically avc/o' = (l 



2Rl 



5 



as in Ref.2^ For i?2 = 0.9 and = 0.8, cryc/o" ^ 1.26 
for H = 0. In our numerical calculation, cryc/cr — 1.24 
at iJ = is a little less than 1.26 because, as we men- 
tioned, rjB is finite. However, when we reduce r]B to 
values smaller than 0.05, our numerical results approach 
this analytic value; for example, for rjB = 0.01 we obtain 
cryc/o' — 1.26. 

A more interesting case is the unitary scattering limit, 
where the vertex corrections are much larger at zero 
field. The parameters Af'' and A^2'' are now, Af'' = 
and 4") = )(/(")) with (4")) = 

where {f3[°^) = J deP{e)p[°\e). For R2 = 0.9 and 
i?3 = 0.8, we obtain 



In 



Po 



Srii 2'KVfVg 
7(0)2 



and 



(0) _ 



z^ln (ipo/zo) 



(42) 



(43) 



where zq = —e + i^t{0) as before. Our numerical results 
are given in Fig. 6. The format is the same as Fig. 5. The 
impurity parameter rju — 10 and b^r = 0,0.1. Referring 
to the bottom frame we see that the application of a 
magnetic field reduces the vertex corrections rapidly as 
the magnetic field increases up to a few times 700. 

For completeness we have also considered the vertex 
corrections to the thermal conductivity of a d-wave su- 
perconductor in the vortex state. It has been shown that 
the vertex corrections are negligible at zero fieldr^ In the 
presence of a magnetic field we obtain 



KVC 



1 + (Vg/Vff 




f deVie)Re 


\ ifie) 




J 


1-Af 


1-4^) 


1 + [Vf/Vgf 




[ deV{€)Re 






J 


[1-4^) 


1-4'^ J 



(44) 



where 



(3), 



z*F{z) - zF{z*) 



and 



[z -f z*){z — z*) 
dF{z) F{z)' 



dz 



(45) 



(46) 



Later we will give A^'-p and A^2^ for the Born and unitary 
limit instead of their general expressions. For isotropic 
scattering, A'^-* and A^2^ vanish again and the thermal 



conductivity becomes 

K _ 1 
Koo 2 

Since /f^(e) = -1 and 



deVU)Re 



r(i)r 



7t(0) 



ItiO) 



arctan 



7t(0) 



(47) 



(48) 



we reproduce the bare bubble result for the thermal con- 
ductivity. 

In the Born limit, we obtain ^j"^-* = a^i\l2^^), A^i'' 

4"'(/f 



„(")/7-(3)\ ^(1) 



-4^\ and A[^'^ 



-Af\ Further 
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FIG. 5: The dc conductivity avc as a function of Eh /'yoo 
including the vertex corrections in the Born limit {jjB ~ 0.05). 
The top frame shows avc/o'oo while the bottom frame gives 
avc i- 6- the ratio of its corrected value (avc) to the dc con- 
ductivity (a) without vertex corrections. The solid (dashed) 
curves in both frames are for — (6^ ~ 0.1). 



simple analysis can show that 1 ± a\ {I2 ) ~ 1 and 
1 ± af\l[^^) « 1 since af^ - 1/ ln[po/7(0)] and e < po 
because of V{e). In the unitary limit, we obtain ^2^-' = 
= and Ai^) = A« = (M^^X/P^ 

Since {(3^°^) -> as po -> 00, 1 - {f3^°^) ~ 1 and 

1 ~ {^2'^) {^2^'') ~ 1- We conclude therefore that the ver- 
tex corrections to the thermal conductivity remain neg- 
ligible even in the presence of a magnetic field; Kyc — 
In Fig. 7 we plot the Lorcnz number, Lyc/Lqq, including 
vertex corrections, vs. the normalized magnetic energy. 
We have used the results for the thermal conductivity 
without vertex corrections since they are small, as ex- 
plained above. While vertex corrections become less im- 
portant at high fields, it is clear form these results that 
a violation from the WF law remains at all fields. 



V. CONCLUSIONS 



FIG. 6: The dc conductivity avc as a function of Eh /'yoo 
including in the unitary limit {jiu — 12). The top frame shows 
o"vc/o"oo while the bottom frame gives avc i-e. the ratio of 
its corrected value {avc) to the dc conductivity [a) without 
vertex corrections. The solid (dashed) curves in both frames 
are for b^, = G {b„ = 0.1). 

Fig.7 (Kim et al) 
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FIG. 7: The WF ratio, with vertex corrections as a function 
of Eh Moo, for the cases considered in Fig. 5 and Fig. 6 with 
bK ~ be- The top (bottom) frame is for the Born (unitary) 
limit. In both cases a violation remains at all magnetic fields. 



We have found that the Wiedemann-Franz law is vio- 
lated in the mixed state of a d-wave superconductor at 
zero temperature at the bare bubble level. Recall that in 
this approximation there is no violation in the absence 
of a magnetic fieldi^ The deviation from the universal 
value of the Lorcnz number Lqq increases with increasing 
magnetic field. When Andreev scattering from vortices 
is neglected, the high field saturated value of the Lorenz 
number can be as low as Loo/2. With Andreev scatter- 
ing, it can be even smaller. We have also found that 
application of a magnetic field affects the effective impu- 
rity scattering rate 7 at zero temperature. In the Born 



limit, 7 increases while in the unitary limit, it decreases 
as is expected. Physically, in both cases, 7 should tend 
toward its normal state value as the field approaches Hc2- 
In the superconducting state with zero field, the relevant 
scattering rate (as is well-known) is much smaller than its 
normal state value in the Born limit while the opposite 
holds in the unitary case. 

We have also considered the vertex corrections arising 
from anisotropy in the impurity scattering. It has been 
shown that in the absence of a magnetic field the vertex 
corrections can be large for the zero frequency electrical 
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conductivity at zero temperature while they are small for 
the thermal conductivity. We found that application of 
the magnetic field reduces the contribution of the ver- 
tex corrections to the electrical conductivity and they 
become unimportant. For the thermal conductivity, the 
corrections remain negligible as in the absence of a mag- 
netic field. 
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